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Abstract 

The cargo motion in living cells transported by two species of motor protein with different intrinsic 
directionality is discussed in this study. Similar to single motor movement, cargo steps forward and 
backward along microtubule stochastically. Recent experiments found that, cargo transportation by two 
motor species has a memory, it does not change its direction as frequently as expected, which means that 
its forward and backward step rates depends on its previous motion trajectory. By assuming cargo has 
only the least memory, i.e. its step direction depends only on the direction of its last step, two cases of 
cargo motion arc detailed analyzed in this study: (I) cargo motion under constant external load; and 
(II) cargo motion in one fixed optical trap. Due to the existence of memory, for the first case, cargo can 
keep moving in the same direction for a long distance. For the second case, the cargo will oscillate in 
the trap. The oscillation period decreases and the oscillation amplitude increases with the motor forward 
step rates, but both of them decrease with the trap stiffness. The most likely location of cargo, where the 
" ^ ■ probability of finding the oscillated cargo is maximum, may be the same as or may be different with the 

trap center, which depends on the step rates of the two motor species. Meanwhile, if motors are robust, 
. i.e. their forward to backward step rate ratios are high, there may be two such most likely locations, 

located on the two sides of the trap center respectively. The probability of finding cargo in given location, 
the probability of cargo in forward/backward motion state, and various mean first passage times of cargo 
to give location or given state are also analyzed. 

Introduction 

Motility is one of the basic properties of living cells, in which cargos, including organelles and vesicles, 
are usually transported by cooperation of various motor proteins [1,2], such as the plus-end directed 
kinesin and minus-directed dynein [3-5]. Experiments found that, using the energy released in ATP 
hydrolysis [6—9] , these motors can move processively along microtubule with step size 8 nm and in hand- 
over-hand manner [10-12]. 

Although numerous experimental and theoretical studies have been done to understand this cargo 
transportation process, so far the mechanism of which is not fully clear. In [13], one basic model is 
presented by assuming cargo is transported by only one motor species and all the motors share the 
external load equally Then in [14], one more realistic tug-of-war model is designed, in which the cargo 
is assumed to be transported by two motor species with opposite intrinsic directionality, and motors can 
reverse their motion direction under large external load. According to some experimental phenomena this 
tug-of-war model seems reasonable [15, 16]. In either of the models given in [13, 14], the only interaction 
among different motors is that, motors from the same species share load equally and motors from different 
species act as load to each other. In [17-19], some complicated models are presented, in which interactions 
among motors are described by linear springs. Recent experiments found that the tug-of-war model might 
not be reasonable enough to explain some experimental phenomena, so several new models are designed 
to try to understand the mechanism of cargo motion by multiple motors [20-26]. Finally, more discussion 
about cargo transportation in cells can be found in [27-35]. 

In recent experiment [36], by measuring cargo dynamics in optical trap, Leidel et at found cargo 
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motion along microtubule has memory. Cargo is more likely to resume motion in the same direction 
rather than the opposite one. This finding implies that, cargo location in the next time depends not only 
on its present location but also on how it reaches the present location. The behavior of cargo depends on 
its motion trajectory, which is different from the assumptions in previous models. In this study, one model 
for cargo motion with memory will be presented. But for simplicity, we assume that the cargo has only a 
little memory, it can only remember the motion direction in its last step. The description and theoretical 
analysis of the model with memory will be first given in the next section, and then corresponding results 
will be presented in the following section. Results will be summarized in the final section. 

Model for cargo motion with memory 

In this study, the cargo is assumed to be tightly bound by two motor species: plus-end (or forward) 
motors and minus-end (or backward) motors. The forward and backward step rates of each plus-end 
motor are u and w, and the forward and backward step rates of each minus-end motor are / and b. 
Obviously u 3> w but b 3> / when the external load is low, since the intrinsic directionalities of the two 
motor species are opposite to each other, and the intrinsic motion direction of plus-end motor is plus-end 
directed (i.e. to the plus-end of microtubule), but the intrinsic motion direction of minus-end motor is 
minus-end directed (i.e. to the minus-end of microtubule). By assuming that all motors from the same 
motor species share the load equally, we only need to discuss the simplest cases in which the cargo is 
transported by only one plus-end motor and one minus-end motor. For example, if there are k plus-end 
motors, the total external load is F c , the forward and backward step rates of one single plus-end motor 
are u c and w c , and the motor step size is l c . Then these k plus-end motors can be effectively replaced by 
one single plus-end motor with load F = F c /k, step rates u = ku c and w = kw c , and step size la = l c /k. 
Since the experiments in [36] showed that, the number of motors moving the cargo is usually the same 
in both directions, this study also assumes the step sizes of the plus-end motor and minus-end motor are 
the same (note, the step size of single plus-end motor kinesin and step size of single minus-end motor 
dynein are the same Iq s» 8 nm [2,9, 12]). 

This study will mainly discuss two special cases: (I) Cargo moves under constant external load. In 
vitro, this constant load may be applied by one feedback optical trap, or In vivo, this constant load may 
be from the viscous environment with invariable drag coefficient. (II) Cargo moves in one fixed optical 
trap, this case is easy to be performed experimentally, and so the corresponding theoretical results are 
easy to be verified. 

Cargo Motion under constant load 

For the sake of convenience, the cargo is said to be in plus-state n + if it reached its present location 
n by one forward step from location n — 1. Similarly, the cargo is said to be in minus-state n~ if its 
previous step is minus-end directed, see Fig. 1(a) for the schematic depiction. In plus-state, the forward 
step rate is higher than backward step rate u > w, but in minus-state the forward step rate is lower than 
backward step rate f < b. So in plus-state, the cargo is more likely to move forward, but in minus-state, 
the cargo will be more likely to move backward. For example, for a cargo in location n, if its previous 
step is plus-end directed, from either plus-state n + — 1 or minus-state n" — 1 to location n, then in the 
next step the cargo will be more likely to move to location n + 1 (plus-state n + + 1), since the cargo is now 
in plus-state n + and its forward step rate u is higher than its backward step rate w. On the contrary, 
if it got to its present location n from location n + 1 (either from plus-state n + + 1 or from minus-state 
n_ + 1), then in the next step the cargo will be more likely to move to location n— 1 (minus-state n~ — 1), 
since the cargo is now in minus-state n~ and its backward step rate b is higher than its forward step rate 
/. This behavior means that the cargo can remember its motion direction of its last step. 
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Let p, p be probabilities of cargo in plus-state and minus-state respectively, then 

dpj dt = fp — wp = —dp/dt. (1) 
Using the normalization condition p + p = 1 , its steady state solution can be obtained as follows 

P = f/(f + w), p = w/(f + w). (2) 
Let U e ff = up + fp, W e ff = wp + bp, then the mean velocity of cargo can be obtained as follows 

V = (U eff - W eff )l = [(u - w)p + (/ - b)p]l = (uf - wb)l /(f + w), (3) 
where Iq is the step size of cargo. The probabilities that cargo steps forward and backward are then 

Ueff = f(u + W) 

(4) 



P+ 



Ueff + W e ff f(u + W )+w(f + b)' 

w(f + b) 



p_=l-p + = 



f(u + w)+w(f + b)- 



Finally, the external load F dependence of rate u, w, f, b can be given by the following Bell approxi- 
mation [37-40], 



u = u e- eoFlo/kBT , w = w e {l -^ Fla ' kBT , f = f e-^ Fl °l k * T , b = b e^-^ Fl °' k ^ 



(5) 



Where eo and t\ are load distribution factors for the plus-end motor and minus-end motor, respectively. 
ks is Boltzmann constant, and T is the absolute temperature. 

Cargo Motion in one fixed optical trap 

This special case is schematically depicted in Fig. 1(b). For convenience, the center of optical trap 
is assumed to be fixed at location 0. For this case, the potential of cargo depends on its location n. 
The potential difference between location n and location n + 1 is AG„ = n[(n + l)Zo] 2 /2 — n(nlo) 2 /2 = 
n(n + 1/2)Zq. Similar as in [19], at location n, the forward and backward step rates u n and w n of cargo 
in plus-state, as well as the step rates /„ and b n of cargo in minus-state, can be obtained as follows, 



-e AG„/fe B T 



; (l-eo)AG n _ 1 /fe B T ) = j- e -e 1 AG n /k B T ^ ^ = foe (l- ei )AG„_i/fc B T^ /g^ 



Where u,w,f,b are cargo step rates when there is no optical trap and any other external load, which 
satisfy u ^> w,b 3> / . For simplicity, this study assumes that eo, e\ are independent of cargo location n. 

Let p n , p n be the probabilities of finding cargo in plus- state n + and minus-state n~ , respectively. One 
can easily show p n , p n are governed by the following equations 

dp,Jdt = U n _xp n _i + fn-lPn-l - K + w n)Pn, (7a) 

dp n /dt = w n+ ip n+ i + b n+ ip n+1 - (/„ + b n )p n . (7b) 

The steady state solution of Eqs. (7a, 7b) are as follows (for details see Sec. A of the supplemental 
materials) 

(fk + b k )u k 



Pn 



'n—X 

n 

.fe=0 



n 

k=n+X 



(uk+x + w k +i)b k 
(uk + w k )b k -x 



(f, 



k-X 



>k-X 



)u k - 



Po, for n > 1, 
po, for ?i < — 1, 



(8a) 
(8b) 
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Pn = -r-Pn = -r- 

On. b n 



n 

,k=0 


n 

,k=n+l 



(fk + b k )u k 
(uk+i + Wk+i)bk 

(u k + w k )bk-i 
(fk-i + 6fe-i)ufc-i 

U 



Po, for n > 1, 



Po, 



for n < — 1, 



Po 



(8c) 

(8d) 
(8c) 



Where po can be obtained by the normalization condition X)n^oo(P™ Pn) = 1- 

The probability of finding cargo in plus-state is p = ^^=-00 P«' an d the probability of finding cargo 
in minus-state is p = ^£^^_ oa The mean locations of cargo in plus-state and in minus-state are 



<n+> 



E 

n— — 00 



npn/p, 



')= E n />™/p> 



respectively. The mean location of cargo is 

+00 



(n) = ^2 n (P™ + Pn) = P( n+ ) + P( n )• 



(9a) 



(10) 



Specially, for the symmetric cases u = b,w = /, i.e. the cargo is transported by two motors with the 
same step rates but different intrinsic directionality, one can verify that p n = p- n and consequently 
p = p, (n~) = — (n + ), (n) = 0. 

The external load dependence of rates u n , w n , f n , b n [see Eq. (6)] means that, for a cargo towed by 
two motors in one fixed optical trap there are two critical values of the cargo location n, 



ksT u 1 

hi - + - - e 



hlf\ 



W 



k B T f 1 



(11) 



where \x] is the smallest integer number which is not less than x, [x\ is the biggest integer number 
which is not bigger than x. The step rates of plus- end motor satisfy u n > w„ for n < n c+ , and u„ < w n 
for n > n c+ . Similarly, the step rates of minus-end motor satisfy b n > f n for n > n c _, and b n < f n 
for n < n c -. The intrinsic directionality of plus- end motor (u 3> w) implies n c+ > 0, and the intrinsic 
directionality of minus-end motor (b 3> /) implies n, c _ < 0. Generally, the critical values n c+ and n c _ 
are different with the mean locations (n + ) and (n~). 

In the following of this section, various mean first passage time (MFPT) problems about the cargo 
motion in fixed optical trap will be discussed. 



Mean first passage time to one of the plus-state 

Let t l n and r l n be MFPTs of cargo from plus-state n + and minus-state n~ to plus-state l + respectively, 
then t l n and T l n satisfy [41,42] 



w n T„_i ~ (u„ + w n )t l n + u n t l n+1 = -1, for n ^ I, 

bnJn-\ ~ (fn + ^n) T ri + = — 1) 



with one boundary condition t\ = 0. 
From Eq. (12a) one can easily get 

Un 



j 



4 4+i , forn^Z. 

W n W n 



(12a) 
(12b) 

(13) 
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Substituting (13) into (12b), one obtains 

1 



u n T Wn A u n A 

'n+1 

w n w n w„ 



(fn + b n ) 



U n+ i + w n +i . 



l _n+]_ t l 



W n +1 



n+1 w n+1 n+2 w n+1 



+ fnt n +l 



l.C. 



^n^n (-^n ~t~ ^n)^n+l ^n^n+2 — ^n 



where 



(«n + Wn)&n 



(/n +f»nK 



-1 



Wn+1 



C^n /n ~t~ ^n 
n — t. 

Wn W n+ i 



(14) 
(15) 

(16) 



Note, Eqs. (14, 15) are established for n ^ I — 1,1. 
Meanwhile, from Eq. (12b) one can get 



fn + K 



n+1 r n p n — 

Jn Jn 

and then by substituting Eq. (17) into Eq. (12a) one obtains 



WnT^i - (ti„ + W„ ) 



/n-1 +b n -l 



fn- 



'n-1 



b n -l 
/n-1 



1 



'n-1 



fn 



fn + b 



fn 



-t! 



fn 



n-1 



1 

fn 



i.e. 



where 



Bn T n-2 — (B n + Fn)T n -i + F n T n — C„, 



®n fan 



fn 



- _ U n U n + W n 

n ~ T ~~f 

Jn Jn—1 



(17) 

-1, 
(18) 

(19) 
(20) 



Eqs. (18, 19) are established for n ^ I. 

The procedure of getting MFPTs t n , t„ is as follows. (1) Getting t l n for n < I - 1 by Eq. (15) and 
boundary condition t\ = (see Sec. B of the supplemental materials). (2) Getting r n for n < I — 2 by 
Eq. (13). (3) Getting t[_ 1 from the special case of Eq. (12b), i.e. &j_it/_ 2 — (fi-i + bi-^rj^ = — 1. 
(4) Getting r n for n > I by Eq. (19) and boundary value t\_ x obtained in (3) (see Sec. C of the 
supplemental materials). (5) Getting v n for n > I + 1 by Eq. (17). This procedure can be summarized 
as follows 



(21) 



n=/-l 



Mean first passage time to one of the minus-state 

Let i l n and f n be the MFPTs of cargo from plus-state n + and minus-state rT to minus-state l~ , 
respectively. Similar as the discussion in Sec. , the MFPTs v n and f n satisfy the following equations 

Wnf n _x - (u n + w n )i n + u n t l n+1 = -1, (22a) 

6nT„_i - (/„ + b n )f l n + fj n+1 = -1, for n ^ I, (22b) 
with one boundary condition f\ = 0. From Eq. (22a) one can easily get 



-l _ u n + W n p U n -i 



y- 1 

'n+1 



1 

it',. 



(23) 



G 



Substituting (23) into (22b), one obtains 

1 



b„ 



i.e. 



^n J n+1 
W n W n W n 



-(/„ + &„) 



-t 



n+l l ri+2 

Wn+1 »n+l 



Wn+1 



with B n ,F n ,C n given by Eq. (16). Note, Eqs. (24, 25) are established for n ^ 
Meanwhile, from Eq. (22b) one can get 



ri fn + b n -i b n 

<Wi = — - — T - - ~ T - 



fn n fn 

and then by substituting Eq. (26) into Eq. (22a) one obtains 



1 

fn 



r , forn^Z, 



fntn+l — 1) 

(24) 
(25) 



(26) 



W n T l n ^ 1 - (Un + W n ) 



fn-1 +b n -l 



fn- 



b n -l 
fn-1 



'n-2 



fn 



fn + b 



™ =l Un =l 

T„ T~T„ 



£ 71 £ 

J n J n 



1 

fi i 



i.e. 



B n f n _ 2 - (B„ + F„)<_! + F n f l n = C n , 



-1, 
(27) 

(28) 



with B n ,F n ,C n given by Eq. (20). Eqs. (27, 28) are established for n ^ I, I + 1. 

The procedure of getting MFPTs t l n , f l n is as follows. (1) Getting f l n for n > I + 1 by Eq. (28) and 
boundary condition f\ = (see Sec. D of the supplemental materials). (2) Getting t l n for n > I + 2 by 
Eq. (26). (3) Getting from the special case of Eq. (22a), i.e. — (w;+i + wi +1 )i l l+1 + ui + ii l l+2 = —1, 
(4) Getting t l n for n < I by Eq. (25) with boundary value t~i +1 obtained in (3) (see Sec. E of the 
supplemental materials). (5) Getting f n for n < I — 1 by Eq. (23). This procedure can be summarized 
as follows 

Eq- (28) , , Eq. (26) rf , ^ , . n Eq. (22a) -j Eq. (25) . - Eq. (23) -j , n 

f/=0 n=Z+l T 

(29) 

Mean first passage time to one given location 

Let 7/ be the MFPT of cargo from state s to location I (either plus-state l + or minus-state l~), then 
one can easily show that 

for s = k^ 
for s = k' 
for s = k^ 
for s = k" 



T = 



f 



t 

fc! 



and k < I, 

and k < I, 

and k > I, 

and k > I. 



(30) 



It is to say that if k < I, a cargo located at k will first reach plus-state l + before reaching minus-state 
l~ . On the contrary, if k > /, it will first reach minus-state l~. Finally, the mean oscillation period T of 
cargo in fixed optical trap can be approximated as follows 



T 



(31) 



see Sec. F of the supplemental materials for its expression. 
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Results 



For cargo motion under no external load, Monte Carlo simulations show that, if the cargo is trans- 
ported by two symmetric motors, i.e., the plus-end motor and the minus-end motor have the same step 
rates, u = b 7 w = f , the cargo will oscillate [Fig. 2(a)]. While for the asymmetric cases, the cargo has 
non-zero mean velocity [see Fig. 2(b)]. On the other hand, if the cargo is put into one fixed optical trap, 
and transported by two symmetric motors, it will oscillate around the trap center with relatively high 
frequency [Fig. 2(c)]. Meanwhile, if the trapped cargo is transported by two asymmetric motors, it will 
also oscillate but its oscillation center may be different with the trap center [Fig. 2(d)]. Both Monte Carlo 
simulations and theoretical calculations show that, for a cargo transported by two symmetric motors and 
put in one optical trap, its oscillation period T decreases with trap stiffness k, motor forward step rates 
u = fe, and motor backward step rates w = f [Fig. 3(a-c)]. Its oscillation amplitude increases with the 
motor forward step rates u = b, but decreases with both the motor backward step rates u = b and the 
trap stiffness k, since high backward step rates and high trap stiffness will prohibit the cargo from moving 
too far from the trap center [Fig. 3(d-f)]. 



Then p is the probability of finding cargo in plus-state, P + is the probability that cargo location n > 
(the center of optical trap is assumed to be at location 0). The meanings of p and P_ are similar. Both 
Monte Carlo simulations and theoretical calculations show that, for a cargo transported by two symmetric 
motors, the ratios p/p and P+/P- are always one, and they do not change with trap stiffness k, forward 
step rates u = b, and backward step rates w = f [Fig. SI]. 

Our results also show that, for cargo motion in optical trap by two asymmetric motors, its oscillation 
period T decreases with trap stiffness k and forward step rate it, but may not change monotonically with 
backward step rate w [Figs. S2(a), S3(a), S4(a)]. But similar as the symmetric cases, cargo oscillation 
amplitude of the asymmetric cases decreases with trap stiffness k and backward step rate w, and increases 
with the forward step rate u [Figs. S2(d), S3(d), S4(d)]. The results in Figs. S3(d), and S4(d) imply 
that, the maximal location n max that cargo might reach toward the plus-end of microtubule depends only 
on the step rates u,w of the plus-end motor, and similarly the minimal location n m ; n that cargo might 
reach towards the minus-end of the microtubule depends only on the step rates b, f of the minus-end 
motor. From the results given in Figs. S2(b,c), S3(b,c), and S4(b,c) one can also see that, different from 
the symmetric cases given in Fig. SI, both the ratio p/p and ratio P+/P- depend on trap stiffness k, 
forward step rate it, and backward step rate w. 

To show more details about the dependence of cargo oscillation on trap stiffness k and motor step 
rates, examples of probabilities p n ,P, and their summation p n + p n are plotted in Fig. 4 and Fig. S5. 
For either symmetric cases or asymmetric cases, the probability profiles are flat for low trap stiffness k, 
indicating that the cargo can reach a farther location from the oscillation center (i.e., with large oscillation 
amplitude) [Fig. S5]. Similar changes can also be found with the increase of motor forward step rates u 
or / [Fig. 4(a, b, d)]. Meanwhile, with the increase of motor backward step rates w or /, the probability 
profile will become more sharp [Fig. 4(c)]. For the asymmetric cases, the most likely location of cargo 
may be different from the trap center [Fig. S5(c)]. One interesting phenomenon displayed in Fig. 4(b, 
d) is that, for either the symmetric cases or the asymmetric cases, when motor forward step rates u, b 
are high, the summation of probability p n + p n may has two local maxima, indicating that cargo motion 
in the positive location (n > 0) is mainly dominated by the plus motor, while its motion in the negative 
location (n < 0) is mainly dominated by the minus motor. 

Let N maxPn , N maxPn , -^(p„+p„) x be the locations at which probabilities p n , p n and their summation 
Pn + Pn reach their maxima, respectively. The results plotted in Fig. 5(a) show that, for symmetric 
motion, N maxfln = — -/V maxpn and their absolute values increase with the forward to backward step rate 



Let 




(32) 



8 



ratio u/w = b/f. The results in Fig. 5(d) show that, for low step rate ratio u/w = b/f, the total 
probability p n + p n has only one maximum which lies at the trap center. However, with increase of these 
ratios, N^ Vn+Pri ^ has one symmetric bifurcation, and its absolute value (see Fig. 4) increases with these 
step ratios. For asymmetric case [see Fig. 5(b)], -/V maxp?i increases with step rate ratio u/w, but iV maxPn 
is independent of it. Which means that, similar as the properties of n max and n m ; n displayed in Figs. 
S3 and S4, -/V maxPii depends only on step rates of the plus-end motor, and N maxPri depends only on step 
rates of the minus-end motor. For asymmetric cases, with the increase of rate ratio u/w, Nr Pn + Pn ) has 
also one bifurcation, see Fig. 5(e). But one of the two values (the negative one) does not change with 
u/w. Which means that, the negative one of ^V(p„+ (9 „) depends only on properties of the minus-end 
motor. Similarly, the positive one of Nr Pn ± Pn ) depends only on properties of the plus-end motor. So 
both the properties of amplitude n max ,n min and the most likely locations N m&xpn , iV maxpn , A r ( Pn+Pn ) 
indicate that, the plus-end directed motion of cargo is mainly determined by the plus-end motor, and 
the minus- end directed motion is mainly determined by the minus- end motor, which is one of the main 
differences with other tug-of-war models [14,18,19,21], and this result is consistent with the experimental 
phenomena [15,16,36]. Finally, the results in Fig. 5(c) show that, the absolute values of N maxPn , N maxPn 
decrease with trap stiffness k, and Fig. 5(f) shows j^( p „+ p „) does not change with stiffness k. So trap 
stiffness can change the oscillation amplitude and the oscillation period (see Figs. 3, S2, and S5), but 
will not change the most likely location Nr Pn+Pn ) of the cargo. Further calculations of probabilities 
p, p show that, for the symmetric cases both p max = p min and (p + p) m i n decrease with step rate ratio 
u/w = b/f, and increase with trap stiffness k [see Figs. S6(a,d)]. Since with large rate ratio u/w = b/f 
and small stiffness k, the cargo will oscillate with large amplitude. For the asymmetric cases, p max ^ p m in, 
Pmax decreases but p min increases with the step rate ratio u/w (i.e. with the increase of the directionality 
of the plus-end motor). Since with large rate ratio u/w, the plus-end motor has high directionality, and 
so the cargo moves fast in the plus-state, which means that the probability p n will be flat with large u/w. 
The plots in Fig. S6(c) show that, although the total probability p n + p n has two maxima, with the 
change of rate ratio u/w, the most likely location of cargo may change from one side of the trap center 
to another side. 

Finally, several examples of MFPTs t l n ,T l n ,i l n ,? l n are plotted in Fig. 6(a,b) and Figs. S7, S8(a,b), 
S9-S12, and examples of MFPTs T l n ± are plotted in Fig. 6(c,d) and Fig. S8(c,d). If m < n < I, then 
ii < t l m ,r l m < T l n , t l n < t l m ,f l m < f < , and T/ l+ < T/ n+ , T l m . < 72- • If I < n < m, then t l n > t l m ,r l m > r l n , 
t l n > *L' T m ^ r ni an d 7^+ > 7^+i XL- > . Moreover, if the trap stiffness k is high and the motor 
step rate ratios u/w and b/f are large, then t l m < T l m ,i l m < f m , T/ n+ < 7^_ for m < n < I, and 
t l m > T [ m ,i\ n > f l m , > T/ n - for I <n< m, sec Fig. 6(a,c,d) and Figs. S7(a,b), S8(c,d),S9, S10(a), 
Sll(b,c,d), S12(a). 

Concluding Remarks 

Recent experimental observations by Leidel et al. [36] show that, in living cells cargo moves along 
microtubule with memory, i.e., its motion direction depends on its previous motion trajectory. In this 
study, such cargo transportation is theoretically studied by assuming that the cargo has the least memory, 
i.e. its motion direction depends only on its behavior in its last step. The cargo will be more likely to step 
forward/backward if it came to its present location by one forward/backward step. Two cases are mainly 
discussed: (I) cargo moves under constant load, and (II) cargo moves in one fixed optical trap. For each 
cases, two kinds of motion are addressed: (i) symmetric motion, in which cargo is transported by two 
species of motor protein which have the same forward/backward step rates but with different intrinsic 
directionality, (ii) asymmetric motion, in which cargo is transported by two species of motor protein with 
different forward/backward step rates. For the symmetric motion (i) of case (I), the mean velocity of 
cargo is zero. But, due to the existence of memory, cargo can move unidirectionally for a large distance 
before switching its direction. One can easily understand that, for the asymmetric motion (ii) of (I), 
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the directionality of cargo with memory is better than that in the usual tug-of-war model by two different 
motor species [f4, 19,21]. For the motion in one fixed optical trap, i.e. case (II), cargo will oscillate. For 
the symmetric motion (i), the oscillation center is the same as the trap center, but for the asymmetric 
motion (ii) , this oscillation center is generally different from the trap center. Usually the oscillation 
period decreases with the trap stiffness k and motor step rates. Meanwhile, the oscillation amplitude 
decreases with trap stiffness k and motor backward step rates w, f, but increases with motor forward step 
rates u, b. The probability p n + p n of finding cargo at location n may have only one maximum, which is 
the same as the trap center for symmetric motion (i) but different with the trap center for asymmetric 
motion (ii). Meanwhile, the probability p n + p„ may also have two maxima. For symmetric motion (i), 
these two maxima are located symmetrically on the two side of the trap center, and their corresponding 
values of probability p n + p n are the same. However, for the asymmetric motion (ii), these two maxima 
are generally not symmetrically located around the trap center, and their corresponding probabilities may 
be greatly different. With the change of ratio of motor forward to backward step rates, the maximum 
with the larger value of probability p n + p n may transfer from one side of the trap center to another 
side. This study will be helpful to understand the high directionality of cargo motion in living cells by 
cooperation of two species of motor protein. Meanwhile, more generalized model can also be employed to 
discuss this cargo transportation process, in which the cargo is assumed to have long memory, its forward 
and backward step rates depend on how long it has kept moving in its present direction. 
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Tables 



Table I. The values of rates u, w, f, b (in unit s _1 ) and optical trap stiffness k (pN/nm) used in the 
plots of Figs. 2-6. The symbol * means that the corresponding parameter is not used in the plot, and 
symbol / means this parameter is one variable in the corresponding plot. Other parameters used in the 
plots arc eo = t\ = 0.5, Iq = 8 nm, and fc^T = 4.12 pN-nm. The stiffness k of the trap used in recent 
experiment of Lcidcl el al. is around 0.02 — 0.09 pN/nm [36]. 
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Figure 1. Schematic depiction of the model discussed in this study to explain the cargo motion with 
memory, (a) is for cargo motion under constant load, and (b) is for cargo motion in one fixed optical 
trap. At any location n, the cargo may be in two different states, plus-state n + and minus-state n~. 
Cargo in plus-state n + means it reaches location n from location n — 1, while cargo in minus-state 
means it is from location n + 1. For a cargo in plus-state n + , its forward and backward step rates are u 
and w respectively. But for a cargo in minus-state n~ , it has different step rates / and b. For the 
constant load cases (a), u > w and b > / mean that, if the cargo is in plus-state n + it will be more likely 
to move forward to location n+ 1. Otherwise, it will be more likely to move backward to location n — 1. 
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Figure 2. Trajectory samples of cargo motion by two motors under constant load (a, b), and in one 
fixed optical trap (c, d). For the symmetric cases (where the step rates of the plus motor are the same 
as the ones of the minus motor, i.e. u = b, w = /), the cargo will oscillate around its initial location (a). 
While for the asymmetric cases, the cargo will have nonzero mean velocity (b). If the cargo is put in 
one fixed optical trap and transported by two symmetric motors, it will oscillate around the trap center 
(c). But for the asymmetric cases, the oscillation center may be different from the trap center. For 
parameter values used in the simulations sec Tab. I. 
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Figure 3. In fixed optical trap, the mean oscillation period T of cargo decreases with trap stiffness n, 
forward rates u = b, and backward rates w = f (in fact, logT decreases almost linearly with log ft, 
logu = log&, and log it; = log/). The oscillation amplitude n max — n m in decreases with stiffness k and 
backward rates w = f, but increases with forward rates u = b. Here n max and n m ; n are the max and 
min locations that cargo can reaches. The circles and squares are obtained by Monte Carlo simulations. 
In (a, b, c), the solid curves are obtained by formulation (31). The solid lines in (d) are obtained by 
n c+ ,n c _ given in Eq. (11), and the solid lines in (e, f) are obtained by n c+ + 3,n c _ — 3, respectively. 
For parameter values see Tab. I. 
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Figure 4. Samples of probability p n and p n for finding cargo in plus-state and minus-state. For the 
symmetric cases probabilities p n and p n are mirror symmetry to each other (a, b, c). Their sum 
Pn + Pn, the probability of finding cargo at location n, might has one maximum [at the center of optical 
trap, see (a, c)] or two symmetric maximum [see (b)]. (d) is one sample for the asymmetric cases. For 
parameter values see Tab. I. 
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Figure 5. The location TV, 



max p 



max p n 1 ^(Pn+Pn ) 



that probabilities p n , p n and their summation 



Pn 



p n reach their maximum. With the increase of rate ratio u/w — b/ f both N n 



and N n 



leave far away from the trap center (a), (b) implies that N maxPn increases with ratio u/w, but N maxPn 
is independent of it. With the increase of trap stiffness k, both N mayiPn and iV maxPn come close the the 
trap center (c). (d, e) show that, with the increase of rate ratio u/w — b/f or rate ratio u/w only, the 
number of maximum of probability p n + p n of finding cargo at location n may change. But (f) implies 
that -^(p„+p„) m , x i s independent of trap stiffness k. For parameter values see Tab. I. 
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Figure 6. Samples of MFPTs t°,r° to plus-state 0+, MFPTs t£,f° to minus-state 0" (a, b), and 
MFPT 7^° from state I to location (c, d). For high trap stiffness k, t„ <0 < r^ l<0 < t ; > < t° >0 for 
MPFTs to plus-state + , and symmetric relations hold for MFPTs to minus-state 0~, see (a). But for 
low trap stiffness, all MFPTs £°,T°,t^,T° increases with the distance between n and trap center 0, see 
(b). Which means that, for different trap stiffness k, the trajectories of cargo from state n + or n~ to 
state + or 0~ are different, (c, d) are MFPTs for one cargo (transported by two asymmetric motors) 
from state n + or n~ to location (plus-state + or CP) and location 1 (plus-state 1 + or 1 _ ). The 
MFPT 7^ is obtained by formulation (30). For parameter values see Tab. I. 



